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Abstract 

The gas-kinetic foundation of fluid-dynamic traffic equations suggested in pre- 
vious papers [Physica A 219, 375 and 391] is further refined by applying the the- 
ory of dense gases and granular materials to the Boltzmann-like traffic model by 
Paveri-Fontana. It is shown that, despite the phenomenologically similar behavior 
of ordinary and granular fluids, the relations for these cannot directly be trans- 
ferred to vehicular traffic. The dissipative and anisotropic interactions of vehicles 
as well as their velocity-dependent space requirements lead to a considerably dif- 
ferent structure of the macroscopic traffic equations, also in comparison with the 
previously suggested traffic flow models. As a consequence, the instability mecha- 
nisms of emergent density waves are different. Crucial assumptions are validated 
by empirical traffic data and essential results are illustrated by figures. 

PACS numbers: 47.50.+d,51.10.-hy,47.55.-t,89.40.-hk 

Key Words: Kinetic gas theory, macroscopic traffic models, traffic instability, dense 
nonuniform gases, granular flow 
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1 Introduction 



Modelling of traffic flow on highways presently attracts a rapidly growing community 
of physicists. Recent research focusses on microsimulation models as well as on fluid- 
dynamic Navier-Stokes-like models p, ^ Q and their derivation from the "microscopic" 
behavior of driver- vehicle units via gas- kinetic equations 0]. Of particular interest is 
the description of traffic instabilities (cf. Fig. ^ which lead to the emergence of stop- 
and-go traffic above a certain critical vehicle density. This phenomenon is illustrated by 
Figure ||. 
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Figure 1: Illustration of the fundamental diagram Q{p) = pVc{p) describing the equilib- 
rium flow-density relationship ( ) and of the temporal course Q{r,t) = p{r,t)V{r,t) of 

the empirical flow ( — ). Obviously, the equilibrium relation Q{r,t) = Qc{p{r,t)) is only 
fulfilled at small densities. After the maximum flow is reached, the temporal development 
of the flow Q(r, t) shows the typical hysteresis phenomenon discovered by Hall |Q which 
indicates a phase transition from almost homogeneous flow to stop-and-go traffic. (Em- 
pirical data: Cross-section r = 41.8 km of the Dutch two-lane highway A9 from Haarlem 
to Amsterdam at November 2, 1994, between 6:30 am and 10:00 am.) 
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Figure 2: Temporal evolution of the mean velocity V{r,t) at subsequent cross-sections 
of the Dutch highway A9 from Haarlem to Amsterdam at October 14, 1994 (five minute 
averages of single vehicle data). The prescribed speed limit is 120km/h. We observe 
a breakdown of velocity during the rush hours between 7:30 am and 9:30 am due to the 
overloading of the highway at r = tq := 41.8km (■■■). At the subsequent cross-sections 

the traffic situation recovers ( : r = ro + lkm; — : r = ro + 2.2km; — : r = ro+4.2km). 

Nevertheless, the amplitudes of the small velocity fluctuations at become larger and 
larger, leading to so-called stop-and-go waves, i.e. to alternating periods of decelerating 
and accelerating traffic. 
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In section |^ it will be shown that most proposed macroscopic traffic models can be viewed 
as special cases of the continuity equation 

dp do 

^ + ^ = with Q = pV (1) 

for the spatial density p{r, t) per lane and and a velocity equation of the form 

dV dV IdV 1 

Here, VdV/dr is the so-called convection term (transport term) which describes velocity 
changes arising from the motion with mean velocity V. The term containing the traffic 
pressure V is the anticipation term and takes into account that driver-vehicle units react 
to the traffic situation in front of them. The relaxation term (\4 — V)/t reflects the 
adaptation of the mean velocity V{r, t) to the density-dependent equilibrium velocity 
(cf. Fig. 1^) with a relaxation time t. 

Apart from the relaxation term, the continuity equation and the velocity equation have 
the form of the fluid-dynamic equations for compressible gases. The additional relaxation 
term reflects that the mean velocity of vehicles decreases at bottlenecks, in contrast to 
ordinary fluids which flow faster. Later on, it will turn out that the relaxation term is 
responsible for the emergence of stop-and-go traffic. 

It should be remarked that all macroscopic traffic models of the above type base on 
phenomenological reasoning in a number of points. Therefore, they imply at least one 
of the following inconsistencies [|12|, |^, ||, ||, depending on the respective model (details 
later) : 

• The variation of individual vehicle velocities is neglected. Some models take it into 
account by the dependence V = pQ of the traffic pressure V on the velocity variance 

e. 

• The velocity variance does not decrease with growing density. 

• The variance does not vanish when the mean velocity vanishes. Therefore, large 
pressure gradients, which can cause the development of negative velocities, are pos- 
sible (cf. Eq. d)). 

• Vehicles are implicitly treated as point-like objects. 

• For a certain density range, the traffic pressure decreases with growing density. This 
is connected with an acceleration of vehicles into regions with larger density (cf. Eq. 
(0)), so that drivers would race into existing traffic jams. 
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Figure 3: Illustration of different suggestions for the equilibrium velocity-density relation 

V^(p) ( — : Helbing 0; : Kerner and Konhauser 0; : Kiihne [|10| and Cremer 

||TT|), compared with empirical traffic data of the Dutch highway A9 from Haarlem to 
Amsterdam, where a speed limit of Vq = 120km/h applies (O: October 14, 1994; +: 
November 2, 1994). The empirical data are mean values of all one minute averages which 
fall into the range [p — 0.5 veh/km, p + 0.5 veh/km). V^(0) < Vq is because of the law that 
speedometers must show the actual speed or more, but never less. 
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• In the course of time, the density p(r, t) can, on certain conditions, exceed the 
bumper-to-bumper density pbb = l/^o, where Iq denotes the average vehicle length. 

• The region of unstable traffic flow, which is related with stop-and-go traffic, is not 
correctly described. In reality, traffic is only stable at very small densities (free flow) 
and extreme densities (slow-moving traffic). 

• Emerging density waves develop a shock-like (i.e. almost discontinuous) structure. 

• The viscosity terms which are implicitly (by the integration method) or explicitly 
introduced for a numerical solution of the macroscopic traffic equations have no 
theoretical foundation. 



In this paper we will show, how to derive a consistent model from a suitable gas-kinetic 
(Boltzmann-like) model. For this purpose, we will extend Paveri-Fontana's Boltzmann- 
like model (cf. Sec. ^ by diffusion effects due to imperfect driving (cf. Sec. Moreover, 
we have to take into account the space requirement by vehicles. This will be done anal- 
ogously to the theory of dense gases and granular flows [|13| (cf. Sec. It will turn out 
that the structure of the resulting fluid-dynamic traffic equations changes considerably 
compared to Eqs. (0) and (0), as well as compared to the equations for dense gases or 
granular materials (cf. Sec. |5.1|) . Reasons for this are the loss of momentum conservation, 
the anisotropy of vehicular interactions, and the velocity-dependence of vehicular space 
requirements. Consequently, a direct transfer of relations for ordinary gases or fluids as 
proposed by previous traffic models is not possible. 



2 Historical evolution of macroscopic traffic models 

From equations (|l|) and (|]) we obtain the traffic models suggested by other researchers, 
if V, Ve, and r are specified in a suitable way: 



In the limit r — > we have 

V{r,t)^VMr,t)). (3) 

The resulting model 

dp / , dVA dp 
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was independently proposed by Lighthill and Whitham and Richards |T^. It 
describes the formation of kinematic waves which propagate with velocity 



c(p) = p- 



dp 



(5) 



relative to 14 . Nevertheless, the model cannot descibe the emergence of stop-and-go 
traffic, since wave amplitudes are not amplified. The wave profile only becomes 
steeper and steeper in the course of time, so that it builds up a shock-like struc- 
ture (i.e. it becomes discontinuous). Consequently, a simulation of Eq. (H) is very 
problematic. Numerically treatable variants of the above model were developed by 
Daganzo . 



Payne 17] suggested the model 



dV ^dV ^ dV 
dt dr dt 



T 



dp/dr dVe 
2p dp 



pr or T 



(6) 



with u = ^\dVc/dp\, corresponding to the traffic pressure 

V = -VJ{2t) . 



(7) 



He derived his model from a microscopic follow-the-leader model [|T8| by means of 
a Taylor expansion [|T^. This model was also used by Papageorgiou [|T^. Cremer 
slightly modified the anticipation term by a factor p/(p + x). Finally, Smulders 



introduced additional fiuctuation terms in the continuity and velocity equation. 
However, the resulting equations still predict the development of shock-like density 
changes, if no "numerical viscosity" is introduced. 

In the limit r — 0, Payne's model implies 

nM) = KWr.<))-I^^. (8) 

p op or 

This results in an additional diffusion term to the Lighthill- Whitham-Richards 
model. 



dp f dVe\ dp _ ^d'^P 
dt \ ^ ^ dp ) dr dr'^ ' 

which solves the shock-formation problem. The diffusion function is 



(9) 



V = T- 



dV 

dp 



dp 



(10) 
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A discrete cell model which is very suitable for real-time simulations of large traffic 
networks was proposed by Hilliges and Weidlich This consists of the discrete 
continuity equation 

^ + ^[Q(M)-Q(^-l,t)]=0. (11) 

Due to the anticipatory driver behavior the flow is modelled by 

Q{t,t) :=p(z,t)K(p(^ + l,t)), (12) 

which assumes that drivers adapt to the velocity in the next cell. Therefore, the 
cell length Ar ^ 100 m must be chosen in agreement with the driver behavior. By 
Taylor approximation of the Hilliges- Weidlich model, we again obtain equation (P), 
but this time the diffusion function is given by 



Ar , 



dV, 



dp 



(13) 



In the special case 



Veip) = Vo[l-^] , (14) 



equation (H) can be transformed to the analytically solvable Burgers equations [E2 



dt ' dr dr 



+ Cir,t)—=V^ (15) 



C{r,t) = Vo{l-^^] . (16) 



with 



Since the stationary solutions of Eq. (^ are stable with respect to fluctuations, the 
above model cannot describe the formation of stop-and-go waves. Therefore, an 
extension by a dynamic velocity equation of the form (^ with V = was proposed 



An alternative approach to the model of Payne was suggested by Phillips [^, who 
derived his model from a Boltzmann-like traffic equation. For the traffic pressure 
he obtained the gas-kinetic relation 

V = pQ, (17) 
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where he assumed 



e(p) = 60 1 



p 



since the velocity variance B should vanish at the maximum traffic density Pmax- 
However, according to this formula the density-gradient dV /dp of the traffic pressure 
will be negative in a certain density range (cf. Fig. p!0|). 

Kiihne as well as Kerner, and Konhauser |]5| avoided this problem by assuming 
0(p) = 00- However, G cannot be interpreted as velocity variance, then, so that 
their equations are not compatible with gas-kinetic traffic models. Moreover, in 
order to smooth out developing shock structures, Kiihne introduced an additional 
viscosity term vd'^Y jdr'^ |Q. In analogy to the Navier-Stokes equations for ordinary 
fluids, Kerner and Konhauser assumed z>(p) = tiq/ p, so that their model corresponds 
to the effective traffic pressure 



V = V,- r,^ — 
or 



with 



= pOo . 



(19) 



Simulations of their model show the formation of density clusters and stop-and-go 
waves at moderate densities The instability region can be determined by a 
linear stability analysis about the stationary and spatially homogeneous solution 
p{r,t) = pe and V{r,t) = Ve(pe)- Inserting the small overall perturbation 



Sp{r,t) := p(r,t) - pc 
SV{r,t) ■.= V{r,t)-VM 



dk p{k) exp[ifcr + (A — iuj)t\ , 
dk V{k) exp[iA;r + (A - iuo)t] 



(20) 



into Eqs. (|T]) and (0), applying a Taylor expansion, neglecting non-linear contribu- 
tions, and applying the orthogonality relations for the complex exponential functions 
yields the following linear eigenvalue problem: 



-A 



ik dV, ^ 1 dV, 



-A 



-ikpe 
r]ok^ 1 



\ Pc dp T dp Pc 
Here, we have introduced the abbreviation 



p{k) 



\ 



(21) 



A := A — ia; with uj := to — kVc{pe) 



(22) 
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k is the wave number, X the growth parameter, and u the oscillation frequency of the 
perturbations. Eq. (^) is fulfilled for the solutions of the characteristic polynomial 

'r]ok^ , 1 



Pe 



This leads to 
A = 



2T ^ 

-±± 

2T 



where we have defined 
1 r]ok'^ 



4T2 



-4^ ± A/3^±i| 
2T ^ ' 



0. 



(23) 



531 



T 



and 



5J : = 



1 



4T2 



Pe 



-a 



1 

+ -, 

r 



1 



op 



dp 



-a 



r dp 

kp, 



9p 



A transition from stability to instability occurs on the condition 

A : 

which implies 



(24) 



(25) 



(26) 



(27) 



(2^ 



Therefore, the equilibrium solution of the Kiihne-Kerner-Konhauser model is unsta- 
ble on the condition 



Pc 



dp 



> 



dp 



1 



rr]ok''' 
pc 



(29) 



This condition is fulfilled at moderate densities, where the equilibrium velocity Ve 
rapidly decreases with growing density (cf. Figs. |^ and ^). According to the con- 
tinuity equation (0), this decrease of velocity causes a further increase of density 
which finally leads to the formation of density clusters. Short wave lengths (i.e. 
large wave numbers) are stable because of the smoothing effect of viscosity t]q. 

The propagation speed of small perturbations relative to V^(pe) is given by the 
relative group velocity 

d 



(30) 
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Figure 4: (a) Illustration of the largest growth parameter A for the instability region (i.e. 
where X{pe,k) > 0). According to the Kiihne-Kerner-Konhauser model (with r = lis, 
rjo = 436km/h, ^/Qo = 54km/h, and the equilbrium velocity Ve{p) depicted in Fig. H) 
traffic flow is unstable at moderate densities and small absolute wavelengths \k\. At 
= we have marginal stability (i.e. A(pe, 0) = 0) due to the conservation of the number 
of vehicles. Instabilities at large absolute wave numbers \k\ (i.e. small wave lengths) 
are surpressed by the smoothing effect of viscosity. Although the model describes the 
instability region qualitatively right, practical experience indicates that stop-and-go traffic 
also develops at much higher densities. 

(b) The representation of the negative group velocity of the unstable mode relative to 
Ve shows that emerging stop-and-go waves move in backward direction, as expected. 
However, the relative propagation speeds seem to be somewhat too large. 



11 



(cf . Fig. ^) . At the transition from stability to instability we find 



c(Pe,A:) = ±iv^=±^^. (31) 

This formula is analogous to the one for the velocity of sound in ordinary compress- 
ible fluids. However, for traffic flows, the plus sign corresponds to the stable mode 
and the minus sign to the unstable mode, so that the forming stop-and-go waves 
move in backward direction. This is in agreement with empirical findings and solves 



a problem raised by Daganzo []T2 



Helbing Q introduced a further dynamic equation for the variance, 
de ^dQ 2VdV IdJ 2,^ 

-Kr + y^ = 7^ — + -ee-e, 32 

ot or p or p or T 

which bases on a gas-kinetic traffic model. This equation is analogous to the 
equation of heat conduction. However, it contains an additional relaxation term 
2(Be — 0)/t. Moreover, despite their theoretical relationship, the variance 6 is 
better not denoted as temperature, and J' is not a heat flow but the flux density of 
velocity variance. 

Another modiflcation suggested by Helbing intended to take into account the flnite 
space requirements s{V) by vehicles. In analogy to the pressure relation of van der 
Waals for a gas of hard spheres (cf. Eq. (|93D), he proposed the relation 

However, in Section it will turn out that this relation cannot simply be trans- 
fered from the theory of ordinary gases to the theory of traffic flow: The dissipative 
and anisotropic interactions of vehicles as well as their velocity-dependent space re- 
quirements change the structure of the fluid-dynamic traffic equations considerably. 



3 Paveri- Font ana's gas-kinetic traffic model 

Prigogine and coworkers EM were the first who proposed the derivation of macroscopic 



traffic equations from a gas-kinetic level of description. However, Paveri-Fontana 
noticed some strange properties of their approach and suggested a modified model which 
will now be discussed. This model describes the temporal evolution of the phase-space 
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density p{r,v,vo,t) of vehicles with desired velocity vq and (actual) velocity v at place r 
and time t, which is governed by the continuity equation 

dp d d f Jv\ d f^dvo\ ( dp\ 



dis 



The acceleration law was specified by 

dv Vq — V 



(35) 



dt T 

delineating an exponential adaptation of the actual velocity v to the individual desired 

velocity vq with a density- dependent relaxation time t{p). Since the individual desired 

velocity vq is usually assumed to be constant, we have 

dvn , , 

— ^ = . 36 
dt ^ ' 

According to Eq. (|3^) , the substantial time derivative of the phase-space density is given 
by temporal changes {dp/dt)^^^ due to discontinuously vnodeWed velocity changes. Usually, 
interaction processes are delineated by this term, since they happen much faster than the 



acceleration processes |^. Paveri-Fontana proposed to use the Boltzmann-like equation 

= (1-p) j dw j dwo\w - v\p{r,v,Wo,t)p{r,w,Vo,t) 

- (l-p) j dw j dwo\v - w\p{r,w,wo,t)p{r,v,vo,t) , (37) 

assuming that a slower vehicle can be immediately overtaken with probability p{p), and 
that the faster vehicle exactly decelerates to the velocity of the slower one, if this is not 
possible. Whereas the first term delineates an increase of p(r, v, vq, t) due to vehicles with 
velocity w > v that must decelerate to velocity v, the second term describes a decrease 
of p(r, V, f 0, t) due to vehicles with velocity v which must decelerate to a velocity w < v. 
According to (pTf), the interaction frequency of vehicles is proportional to their relative 
velocity |f — w| and to the phase-space densities p of the interacting vehicles. Their desired 
velocities vq and wq are not changed by interactions. 

The gas-kinetic equation allows a derivation of macroscopic traffic equations for the spatial 
vehicle density 



p{r,t) ■= j dv J dvop{r,v,vo,t) (38) 

and the velocity moments 

{v'M) ■■= jdvj ^^o^^^t^o)'^^^^!^ , (39) 
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in particular the mean velocity 

V{r,t) := (v) = I dvvP{v;r,t) (40) 

and the velocity variance 

e{r,t) := {{v -Vf) = I dv{v-VfP{v]r,t), (41) 



where we have introduced the velocity distribution 



(42) 



The macroscopic equations are obtained by multiplication of the gas-kinetic equation (|34|) 
with ip{v) := 1, V, v"^ and subsequent integration over v and vq. Finally, one finds the 
fluid-dynamic traffic equations 

dt '''' dr ^ dr ^ ^ ^ 

dV dV IdV 1 

dt dr p dr ~^ T ' ^ ^ 

^ + V— - - — + -Q) (45) 

dt dr p dr p dr T ^ 



which was shown in Refs. We recognize that the density equation (^3|) and the 

velocity equation (^) have again the form of Eqs. (|I|) and (Q). However, we have addi- 



tionally found theoretical relations for the traffic pressure 

V ■= p{{v - Vf) = p& (46) 

and the equilibrium velocity 

Ve:=Vo-T{p)[l-pip)]pe. (47) 



Moreover, we have obtained the variance equation ( ^Sf ) with the flux density of velocity 
variance 

J:=p{{v-Vf), (48) 

the equilibrium variance 

Q.:=C-\r{p)[l-p{p)]p{{v-Vf), (49) 

and the covariance 

C:= {^[v-V){vo-Vo)). (50) 



14 



For the evaluation of J7'(r, t) we need a mathematical expression for the velocity distribu- 



tion P{v; r,t). This can be approximately obtained by Grad's method of moments ^7 
which uses the expansion 

N 



P{v;r,t) ^ ^a„(r,t) 



Qn 



1 



n=0 



exp 



vy 



26 



(51) 



The functions an{r,t) can be expressed by the velocity moments {v'^) {0 < k < N) which 
are governed by the macroscopic equations. In gas theory the macroscopic equations are 
closed after the variance equation (A^ = 2), since higher velocity moments are varying 
on much faster time scales, so that their adiabatic elimination is justified. For vehicular 
traffic as for granular materials, a separation of time scales is not this simple. Therefore, 
one usually restricts to the number of macroscopic equations, which allow to delin- 
eate the instabilities under consideration. Sela and Goldhirsch ||2^ have shown that this 
method yields reliable and "universal" results in the sense that the inclusion of additional 
macroscopic equations gives no fundamental but only minor corrections. 

From our discussion in Sec. ^ we know that the continuity and velocity equation are 
sufficient for a description of emerging stop-and-go traffic. Therefore, we choose = 1, 
so that we must express 6(r, t) in dependence of p{r,t) and V{r,t): 



Q{r,t) := ee(p(r,t),\/(r,t)). 

Because of the conditions (1) = 1 and (v) = V{r,t), we find ao{r,t) = 1 and ai(r, t) 
so that we get the Gaussian distribution 



1 



exp 



vy 



26, 



This leads to 



i7 = and 6e = Ce 



(52) 
= 0, 

(53) 
(54) 



where Ce is the equilibrium covariance (cf. Fig. |^). In Ref. it has been shown that Ce 
is given by the implicit equation 



Ce = ((^^0 - Vof) - 2r(p)[l -p(p)]pCc 




(55) 



Therefore, 9e is only a function of p and not of V (cf. Fig. Eqs. (|5^) and (^) together 
with (|47|) imply the relation 



•Mm 



TC 



(56) 
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In homogeneous traffic situations, the variance of desired velocities Qq := {{vq — Vq)^) 
should be a constant and, in particular, independent of density. Therefore, relation (|56|) 
predicts that the variance is still finite when the equilibrium velocity Vg vanishes. This 
is, of course, a paradoxial result. 
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Figure 5: Comparison of empirical velocity distributions at different densities ( — ) with 
frequency polygons of grouped Gaussian velocity distributions with the same mean value 
and variance ( — ). The class interval lengths are 5 km/h. A signigicant deviation of the 
empirical relations from the respective discrete Gaussian approximations is only found 
at density p = 40 vehicles/km, where the two minute averages of the single vehicle data 
may have been too long due to rapid stop-and-go waves (cf. the mysterious "knee" at 
p ^ 40 veh/km in Fig. |^. The velocity distributions keep their unimodal form even at 
high densities. (Data: Dutch highway A9 with constant speed limit 120 km/h.) 



Another problem was recognized by Shvetsov |2g]: According to Paveri-Fontana's equa- 
tion, vehicles are not decelerated to velocities less than the minimum desired velocity of 
all drivers. Therefore, it cannot describe the development of so-called phantom traffic 
jams, where all drivers want to drive fast, but produce a slowly moving traffic jam. 

As a consequence, we must modify Paveri-Fontana's equation somewhat. This is done in 
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Figure 6: Empirical skewness of the velocity distribution on the basis of mean values of 
one minute averages (O: October 14, 1994; +: November 2, 1994). The absolute value 
of the skewness is rather small and mostly lies between and 0.5. (The variation of the 
data at higher densities results from the small amount of one minute data that could be 
averaged.) 
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the next two sections. 
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Figure 7: Illustration of the empirical equilibrium variance-density relation 0e(p), evalu- 
ated like in Fig. | (O: October 14, 1994; +: November 2, 1994; — : fit function). 



4 Imperfect driving 



Up to now, we have assumed perfectly driving vehicles that are able to keep a constant 
speed and to exactly adapt to the velocity of a slower vehicle in front. Therefore, we will 
replace Paveri-Fontana's interaction term by the more general formula 



9p\ 



dp 
dt 



where 



dp 
di 



diS 



1 



diff 



dp 
dt 



(Vp) 



int 



(57) 



(5^ 



describes velocity fluctuations due to imperfect speed control. The diffusion function V 
must be chosen in such a way that the velocity-diffusion term {dp/dt)dis cannot produce 
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negative velocities. In the following, we will assume 



Div) 



(59) 



in order to obtain a dimensionless relation a{p). According to Eq. (^91), velocity fluctua- 
tions are related to the time scale r of velocity adaptation. Moreover, comparing (^) with 
the formula (|Axp) = VAt for spatial diffusion with diffusion constant V, the absolute 
velocity displacement |Af| is approximately proportional to v during small time intervals 
At, which sounds plausible. 

Now, we will generalize the interaction term. By 



dp 
di 



int 



dw J dv' j d'^' J dwo\w — v'\a{w' ,v\w,v')p{r,w,wo,t)p{r,v' ,vo,t) 

v'>w w'>v 

dw J dv' J dw' J dwo\v — w\a{w' ,v'\w,v)p{r,w,wo,t)p{r,v,vo,t) (60) 

w<v w'>v' 



we can take into account vehicles that do not exactly adapt to the velocity of a slower car 
in front. a{w' ,v'\w,v) denotes the differential cross section and describes the probability 
with which two vehicles have velocities v' and w' after their interaction, if they had 
velocities v and w before. In the following, we will choose 



a(w', v'\w, v) 



P 



I3w 



exp 



w — V 
I3w 



5{w' — w) 



(61) 



where 5{.) denotes Dirac's delta function. According to this formula, the velocity w of the 
slower vehicle is not influenced by the velocity of a faster vehicle with velocity v behind 
it (i.e. w' = w). The exponential function reflects that the faster vehicle slows down to 
a velocity v' < w. Although the velocity v' = w is most likely, smaller velocities v' < w 
also occur with a certain probability due to an imperfect adaptation of velocity. For the 
degree (3 of imperfection we have < /5 <^ 1. In the limit /5 = 0, formula (|6l|) becomes 
a{w' ,v'\w,v) = {l—p)6{v' — w)6{w' — w), so that we get back Paveri-Fontana's interaction 
term (g^)- 

In order to obtain the macroscopic traffic equations corresponding to the generalized 
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gas-kinetic equation, we have to evaluate the interaction terms 



dv I dvoi/j{v) J 

/ int 

p^(r, t) J dv J dw J dv' J dw' ^{v)\w — v'\a(w',v\w,v')P{w]r,t)P{v'] r,t) 

v'>w w'>v 

p^{r,t) dv dw / dv' / dw' 'ijj{v)\v — w\a{w',v'\w,v)P{w;r,t)P{v;r,t) 



p {r,t) J d^ J J J dw' \v — w\a{w' ,v'\w,v) 

w<v w'>v' 

X [^{v') - ^{v)]P{w] r, i)P{v] r, i) , 



(62) 



where we have interchanged variables to obtain the final result. With (|53D and (^) we 
find J(l) = and 



X{v) = -(l-p)p' 



l{{v - Vf 



'l-p)p'\l3VQ, + 2l3'V 




2ir2. / 



(63) 



Since the diffusion term {dp/dt)dm yields no contribution to the density and velocity 
equation, but the contribution 2a{V'^ + 6c)/t to the variance equation, we arrive at the 
corrected equilibrium relations 



V,{p,V) = Vo-r{l-p)p 



e,{p,V) = C, + a{V^ + Q,) + ^{l-p)p{(3Ve, + 2(3^V\l-' 



(64) 
(65) 



Note that the equilibrium relations Ve and 6c not only depend on the density p, but also 
on the mean velocity V, now. The previous formulas result for a = j3 = 0. 



4.1 Determination of the parameters from empirical data 

The calibration of the functions a{p), f3{p), p{p), and r(p) is a difficult problem, since we 
only have empirical relations for V^(p) and 6c (p) (cf. Figs. § and|^). Therefore, we set 
/? = in order to reduce the number of parameters. (Setting a = does not allow to 
describe the finite equilibrium variance at small densities if Cg = 0, see below.) Since the 
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equilibrium covariance Ce is not anymore given by relation (0), we calculate it via Eq. 

Ce(p) = (1 - a)ee(p, V;(p)) - a[V,{p)f . (66) 
This expression is completely determined by K(p) and 6e(p), if we define 

9e(p,K(p)) , . 

p-Zx[K(p)]2 + ee(p,K(p))' ^ ^ 

Consequently, the equilibrium variance 

e.(p,n = ^°'f + (68) 

1 — a 

is given by the diffusion effect q;V^^/(1 — a) at high densities, whereas it is dominated by 
the covariance at low densities, which is very plausible. In the case of speed limits, 
however, we have Ce ~ 0, since all vehicles have approximately the same desired velocity 
fo ~ Vq. In this situation we define 

( \ Qc(p,K(p)) 



The product r(p)[l — p(p)] can be easily obtained via Eq. ( |103| ), which is a generalization 
of Eq. (H): 

PX0c(p,K(p)) 

A determination of the factors p{p) and r(p) is only possible by means of additional 
assumptions. A detailed analysis shows that the relaxation time is given by the relation 

.(P) = ^. (71) 

where g(p) is the proportion of freely moving vehicles and tq ~ 8s the relaxation time of 
vehicles which are not impeded during their acceleration 0. In addition, a complicated 
relation of the form p(p, 5'(p), K(p)) can be derived (cf. Ref. 0). The resulting density- 
dependent functions g(p) = ro/r(p) and [1 — p{p)] are depicted in Figure |. It is very 
important that, although p{p) approaches the value in the limit p Pmax, the relaxation 
time r(p) remains finite, since the maximum density Pmax is somewhat less than the inverve 
of the average vehicle length Iq. 
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Figure 8: Illustration of different relations estimated from empirical traffic data: The 
density-dependent probability (f — p) that a fast car cannot immediately overtake a slower 
one ( — ) starts at zero and reaches the value one in the limit p Pmax- The proportion 
q{p) of freely moving vehicles q{p) is about one for very small densities and approaches a 
finite value with growing density ( — ). Symbols represent the proportion of vehicles that 
have a time headway greater than 2.5T, where T = 0.8 s corresponds to the safe time 
headway (O: October 14, 1994; +: November 2, 1994). 
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5 Dynamics at high densities 



Up to now, vehicles have been imphcitly treated as point-hke objects. Therefore, the 
model must now be corrected for the finite space requirements of vehicles. This will lead 
to further modifications of the interaction term. In order to illustrate the applied method, 
we will first discuss the gas-kinetic and fluid-dynamic description of granular material like 
sand, powder, or pills |jl3|, E9|, BO . 



5.1 Granular flow 

Granular materials have found a broad interest due to their various instability phenomena 
like density waves, avalanches, cluster or heap formation, convection, or size segregation 
In order to work out the similarities with traffic flow, we will focus on the 



description of density waves in sand falling through a narrow pipe [33|. Similar to the 



discussion by Riethmiiller et al. we assume a vertical tube of diameter s, in which 
spherical grains (or cylinders of height s) with radius s/2 are falling. Since grains do 
not have a desired velocity Vq, we are confronted with the phase-space density p{r,v,t) of 
grains with velocity v, this time. The corresponding gas-kinetic equation is 

It exactly corresponds to the gas-kinetic traffic equation in the case of a speed limit Vq, 
where the desired velocities of all drivers approximately agree {vq ~ Vq), so that vq is 
not a variable anymore, but a fixed parameter. Since velocity changes of grains result 
from acceleration due to gravity g as well as from sliding friction at the wall (and by 
displacement of air) with a friction coefficient 7, we have the acceleration law 

^ = g-jv. (73) 

The diffusion term ^Vd'^ p / dv"^ with diffusion constant T> := 2769 describes variations of 
the grain velocities v due to fiuctuating infiuences of the wall (and the displaced air) in 
accordance with the fluctuation- dissipation theorem. 

Multiplying the gas-kinetic equation ( |72D with ip{v) := 1, v, or f^, and integrating over v 
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gives the macroscopic equations 

dt'^ dr ^ dr ^ ^ ^ 

— + V— - --— + { - V) + ^^ (75) 
dt dr p dr p ' 

dt dr p dr p dr p 

where we have taken into account the conservation of the number of vehicles (X(l) = 0). 
A comparison with the fluid-dynamic traffic equations shows that the equations agree 
with each other if we substitute g ^ Vq/t, 7 ^ 1/r, ©o ^ aiV^ + 9), and specify 
the interaction term like in Eq. (|37D with vq = Vq = wo- Although this specification was 
suggested by Riethmiiller et al. |Q , we will apply the theoretical relations for dense gases 
and granular materials [0, instead, in order to include effects due to high densities and 
momentum conservation. The interaction terms have then the form 

T{tp) = J '^^ J J J dw' \v — w\a{w' ,v'\w,v)['iIj{v') ~ 4){v)]/52{r + s,w;r,v;t) 

w<v w'>v' 
w>v w'<v' 

(77) 



In contrast to the formula (^) for vehicular interactions, we have an additional contri- 
bution due to backward interactions of grains (last term). Moreover, we have taken into 
account that pushed grains are located at places r + s ot r — s. Finally, the pair distribution 
function p2 describes velocity-correlations between interacting grains. 

The expression (|7^ can be simplified by interchanging variables: 



^('^) = J J J J 1^ ^ ^'1^' ^) 

x{[i>{v') - i>{v)]p2{r + s,w;r,v;t) 

+ [iIj{w') - ilj{w)]p2{r, w;r - s,v;t)} . (78) 
Applying the Taylor approximation 

d „ 

P2(r + s, w; r, v; t) ^ p2(r, w;r - s, v; t) + s—P^ir + s, w; r, v; t) (79) 

dr 



yields 



I(i.)^Ui.)-?IM (80) 
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with the source-like contribution 



2s(V') = J J J J I"" ~ "^l""^' 

w<v w'>v' 

x{['ijj{v') + ^{w')] - [ip{v) + tjj{w)]}p2{r,w;r - s,v;t) (81) 
and the flux-like contribution 

Xf(ip) = —s J dv J dw J dv' J dw' \v — w\a{w',v'\w,v)['ip{v') — 'ip{v)]p2{r -\- s,w;r,v;t) 

w<v w'>v' 

(82) 

describing colhsional transfer. In this representation we immediately see that the source- 
hke contributions vanish for collisional invariants ip, for which ilj{v')-\-ip{w') = ip{v)-\-il!{w) 
holds. As a consequence, we have only flux-like contributions to the velocity and variance 
equation for ordinary gases or fluids. However, in granular material the conservation of 
kinetic energy gets lost due to inelastic, dissipative interactions. Therefore, a source-like 
contribution is expected in the equation for the so-called granular temperature 6 (cf. Eq. 
(PPP). Since momentum is conserved during granular collisions, these yield only a fiux-like 
contribution to the velocity equation, which results in a corrected pressure relation (cf. 
Eq. (|93|)). Therefore, we have again a velocity equation of the form (^ with a density- 
mdependent equilibrium velocity K = g/l- Consequently, the instability condition for 
the granular flow is given by Eq. (^), if we assume the validity of approximation (|^. 
Due to dVc/dp = 0, the stationary and homogeneous solution of the granular density and 
velocity equation is stable at all densities pe- For this reason, the approximation (|52D is 
not valid for granular flows, and we must apply Grad's method with N > 2. For N = 2 
we obtain the Gaussian-shaped Maxwell-Boltzmann distribution 



P{v; r, t) ^ 1. — exp 
V2^ 
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^ '83) 



because (1) = 1, {v) = V, and {{v — V)^) = 9 imply ao{r,t) = 1 and ai(r, t) = a2(r, t) = 0. 

The explicit calculation of the interaction terms now calls for a specification of the 
differential cross section (t{w' ,v'\w,v) of granular collisions. Since both interacting grains 
have to be equivalently treated (isotropy condition), the collision law v' = fiv-\-fi'w implies 
w' = piw + pi'v. Due to momentum conservation {v' + w' = v-\-w) we must set yu' = (l — p,). 
Therefore, we find the relation 

a{w', v'\w, v) = 5{v' - [fiv + (1 - fi)w])5{w' - [nw + (1 - /i)f]) , (84) 
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which is invariant with respect to interchanging the particles, is a characteristic param- 
eter of the granular material. It is related to the amount of energy dissipation because 
of 



{v" + w'^) = K + ^ ) - 2yu(l - P)iv - ■ (85) 

Since collisional energy must not be produced during the interactions, we find < /i < 1. 
Moreover, since we have v > w before the collision and v' < w' after the collision, we 
obtain the further restriction /i < 1/2. For /i = we have completely elastic collisions 
with an interchange of particle velocities, so that the interaction terms in the macroscopic 
equations ( [75D and (|76|) vanish in the limit s <^ 1/p (like for ordinary gases or fluids) . 
However, due to its analogy with vehicular traffic, we will focus on the extremely inelastic 
case yU = 1/2 in which both particles have the same velocity after their collision. 

Next, we will specify the pair distribution function p2 of interacting grains. This is usually 



expressed by their one-particle phase-space densities p in the following way |T3| : 

P2(r + s, w; r, v; t) = x{r + s/2, t)p{r + s, w, t)p{r, v, t) . (86) 

The factor 



reffects the increase of the particle interaction rate ||3^ because the grains are extended 
by an amount s/2 around their centers, so that they collide earlier. A first-order Taylor 
approximation of p2 gives 

X(r ± s/2, t)p{r ± s, w, t)p(r, f , t) ^ xij, t)p{r, w, t)p{r, f , i) 



x<l±^^±s 
2x or 



Idp v-VdV 1 f{v-Vy \ dQ 



pdr^ Q dr ^ 2Q\ dr 



With this, evaluating the interaction terms X(?/') and neglecting higher order derivatives 
as well as products of derivatives (i.e. neglecting Navier-Stokes and Burnett corrections) 
finally yields the continuity equation and the following Euler-like equations for granular 
fiows: 

dV dV Id 

7^ + ^TT = — + (1 + 3^)Pco„]7^ + 27(00 - 6) - ^(1 - p)^f^ ■ (90) 
dt dr p dr 1 — ps 
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That is, in the velocity equation ( pOf ) the interaction term yields an additional contribution 

ps 



'l-p)V, 



(91) 



1 — ps 

to the pressure. This diverges in the limit p Pmax = 1/s of extreme densities, whereas 
V = p0 vanishes, since the equilibrium variance Qe{p) given by the implicit equation 

4 pQ^ip) 



ee(p) = eo-p(i-/i)- 



(92) 



7v/7r 1 - ps 

vanishes. In the elastic case p = we obtain the formula of van der Waals for the total 
pressure Vtot in a gas of hard spheres: 



tot 



V + Pcorr = P 1 + (1 - p) 



ps 



1 



ps 



ti=0 



V 



1 — ps 



(93) 



The last term in the variance equation ( pOD results from energy dissipation during granular 
collisions. Moreover, we have discovered a new contribution 3pPcorr to the pressure, which 
(to the knowledge of the author) has not been reported elsewhere, presumably because it 
only plays a role in very inelastic cases (p 9^ 0). 

A linear stability analysis about the stationary and spatially homogeneous solution of 
Eqs. ([7^), (p9|), and (pO| ) leads to the characteristic polynomial 

'99, 



2 



- 1 



dVi 



tot 



-2-fk 



dp 
dp 



+ 



1 _ ^® 



(^ot + SpVo 



dQ 



de 



+ 



dVtot dQe 
dQ dp 



0. 



(94) 



Its numerical investigation shows that the Euler-like equations for granular flows are un- 
stable if the noise level Go is sufficiently large (cf. Fig. The mechanism of the clustering 
instability originates from the increase of granular pressure with growing density. This 
causes a reduction of velocity (cf. Eq. (p9D) which results in a further compression (cf. 
Eq. (^). Since Eq. (|9|) implies 

\{pe, k) = \{pe, -k) and uj{pe, k) = u^p^, -k) , (95) 

the relative propagation speed 



u{pe, k) := 
and the relative group velocity 
cipck) : 



pjjpe, k) _ uj{pc, -fc) 
k ~ -k 

duj{pc, k) 



dk 



duj{p^, -k) 
d{-k) 



-u{pc, -k) 



-c{pc, -k) 



(96) 



(97) 



of perturbations with respect to \4 are antisymmetric functions in k. This is in contrast 
to the situation for vehicular traffic (cf. Fig. [1^). 
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Figure 9: Illustration of the largest growth parameter A in units of 7 for the instabihty 
region (i.e. where X{pe,k)/^ > 0). The wave numbers in units of 1/s are restricted to 
the range — 27r < ks < 27r, since wave lengths shorter than the diameter s of the grains 
cannot be propagated by the granular material. Granular flow is stable for wave numbers 
k ^ 0. Emergent density waves appear above a density-dependent critical value of \k\. 
These develop faster with growing absolute wave numbers |A;| and increasing densities pe- 
A smaller degree of dissipation /i (here: = 1/2) results in an expansion of the instability 
region towards higher densities pe and smaller absolute wave numbers \k\. Decreasing 
the noise level ©0 leads to a reduction of the instability region, until it vanishes below a 
certain critical value. 
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5.2 Traffic flow 



Our previously discussed traffic model is now very easily extended to the description of 
dense traffic. Again, we can apply formula ([78]). However, inserting ( ]HT| ) shows that 
the last term, which corresponds to backward interactions, vanishes. Therefore, the only 
difference between (|62[) and (|78|) is the replacement of p'^P{v;r,t)P{w;r,t) by the pair 
distribution function p2{r + s,w;r,v;t). In addition, we must take into account that 
vehicles require a velocity-dependeiat space of about 

s{V):=^ + TV, (98) 

Pmax 

since drivers keep on average a save distance of length TV, where T ^ 0.8 s is about the 
reaction time. As a consequence, we get the modified relation 



P2{r + s, w; r, v; t) = x(r + TV, t)p{r + s{V),t)P{w; r + s{V),t)p{r, t)P{v; r, t) (99) 



with 

r + TV{r, t) is the interaction point of a vehicle at place r with a vehicle at place r + 
s(y{r,t)). In contrast to our discussion of grains, r and r + s{V) are not the centers but 
the fronts of the vehicles. 

Now, we can evaluate the collision terms (ffSl). Applying Taylor approximations for x(r + 
TV, t) and p(r + s{y),w, t), we finally obtain the velocity equation 

dV ^BV dp dV d& , d^p , d^V , d'^Q 1,,, , , 

ot or or or or or^ or-'' or''' t 

We have neglected products of partial derivatives (Burnett corrections), but taken into 
account second order derivatives (Navier-Stokes terms). The abbreviations used in the 
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above equation are: 



ai 



e 
p 



+ (i-p)p 



P 



dp 



e 1 



/3 



+ /3V) 



aa = (l-p)p|sx 
as = 1 + (1 - p)psx 
bi = {l-p)p 



Ti 2 



-TV 



dx 
dV 



6 1- 



f3 



2 4^716 



2p 



2 dp 



(3 



e 



6 1-^ +/?^\/- 



TC 



2 dV 



&3 ^ 

and 



4V7re 



Ve = Vo-{l-p)TpX 



e 1 



/5 



(102) 



(103) 



Again, we will eliminate the variance equation by means of the approximation 0(r, t) ^ 
Gc(p(r, t), y(r, t)). In the case Ce ~ related to a speed limit, the equilibrium variance 
6e(p, V) is given by the implicit equation 



ee(p, V) = a{V^ + Oe) + J(l - p)px + 2/5V^/ ^ 



(104) 



which is obtained by evaluation of the variance equation. Consequently, ©e vanishes when 
V becomes zero, which is required for consistency 0, ^. Applying the above approxima- 
tion for the variance and again neglecting products of partial derivatives (i.e. Burnett 
corrections), we find 



dV ^dV 
dt dr 



p dp dr dr dr"^ p dr"^ r ^ 



(105) 



where 



dV, 



tot 



dp 



p ai + as 



P[h: 



de, 

dp 

dQe 
dV 



de. 



a2 — a-i 



bi + k 



dV ' 

g©c 

dp 



(106) 
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and Ve are functions of p and V. 

It is obvious that the loss of momentum conservation and the velocity-dependence of the 
vehicular space requirements have drastically changed the structure of the velocity equa- 
tion: The terms containing a and b are completely new compared to the Navier-Stokes 
equation for ordinary or granular fluids and compared to all previous macroscopic 
traffic models (cf. Sec. 0). Fortunately, the model solves the problems which were men- 
tioned in Section |l|: 

• The variation of individual vehicle velocities is taken into account by the variance 

• The equilibrium velocity Ve and the equilibrium variance Ge are monotonously de- 
creasing functions that vanish at the finite density Pmax- 

• The equilibrium variance 0e(p, V) vanishes when V{r,t) = 0. 

• In equilibrium (i.e. for p{r,t) = p^ and V{r,t) = Vc{pc)) the density-gradient dV /dp 
of the total traffic pressure V is non-negative, so that the latter is a monotonously 



^max 



increasing function of pe (cf. Fig. |T0|). The singularity of the pressure at p = 
guarantees that the upper bound pmax < phh cannot be exceeded by the density 
p(r,t). 

We were able to explain viscosity as an effect of the finite distance s{y) between 
interacting vehicles. Moreover, we obtained a theoretical expression for the viscosity 
//(p, \^), which not only depends on the density, but also on the mean velocity. 
In equilibrium p{r,t) = pe, V{r,t) = V^(pc), the viscosity is non-negative. The 
singularity of viscosity at p = Pmax causes that extreme changes of V{r, t) and 
p{r,t) are smoothed out, so that the shock-formation problem is solved. 



The fiuid-dynamic traffic model (^3|), ( p.05| ) is able to describe the emergence of 
stop-and-go traffic at medium densities (cf. Fig. |I^a). A linear stability analysis 
leads to the characteristic polynomial 

X^ + X[^ + --ika- -—f + ikpe + -— ^ + = . 107 

\ Pe T T dV J \ Pe dp T dp J 

The transition from stability to instability occurs on the condition 
hpeTk"^ + ( 1 + 



Pe 



dVe 



dp 



dVe 



dV 



rjTk'^ \ ( a , ja^ dVtot 
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Figure 10: Comparison of the density-gradients of the total traffic pressure Vtot ( — ), the 

ideahzed pressure V = p0 of point-hke objects ( — ), and the correction term Vcon ( )i 

estimated from empirical data. For p ^ 20veh/km and large densities, the idealized traffic 
pressure V of point-like objects decreases with growing density. However, this is more 
than compensated by the pressure correction Vcorr due to the finite space requirements of 
vehicles. In particular, the increase of pressure with growing density diverges in the limit 

P ^ Pmax- 
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Figure 11: Illustration of the viscosity function estimated from empirical data. Obviously, 
the viscosity is non-negative, as expected, and diverges at large densities. Viscosity causes 
the smoothing effect which is necessary to avoid shock-like structures and to facilitate 
numerical simulations of the macroscopic traffic equations. 
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In contrast to Eq. (p9D, an additional instability appears at large absolute wave 
numbers \k\ > 130/km. However, these are connected with wavelengths that are 
smaller than the average vehicle distance 1/p so that they cannot be propagated by 
the discontinuous vehicular fluid. Perturbations which can actually be propagated 
by the vehicles are moving in backward direction with respect to Ve (cf. Fig. |T2|b). 

6 Summary and Outlook 

It was shown that almost all macroscopic traffic models can be viewed as special cases of 
the continuity equation and a certain velocity equation. Although these equations have a 
close similarity to the hydrodynamic equations for ordinary fluids, they did not fulffil all 
consistency requirements for realistic traffic flow models. 

Finally, it turned out that the correct structure of fluid-dynamic traffic equations looks 
considerably different. We obtained this result from a refined version of Paveri-Fontana's 
gas-kinetic traffic equation which was extended by the effects of imperfect driving and 
vehicular space requirements. In particular, we have applied the gas-kinetic theory of 
dense gases and granular materials. 

Nevertheless, despite the phenomenologically similar behavior of traffic flow and granular 
material falling through a narrow vertical pipe, the governing equations and instability 
mechanisms are completely different. This was illustrated by numerical results of linear 
instability analyses and originates from the fact that momentum is conserved by granular 
collisions but not by vehicular interactions. 

The parameters and relations occuring in the final model have been estimated from em- 
pirical traffic data. In addition, we have empirically tested the approximations that we 
made during the derivation of the fluid-dynamic traffic model. Finally, it was shown that 
all consistency criteria are met. 

Future work will compare simulation results of the various discussed traffic models. More- 
over, the refined traffic flow model (which treats the highway lanes in an overall manner) 
will be extended to a model for the different interacting highway lanes. 
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Figure 12: (a) Illustration of the largest growth parameter A for the instability region 
of traffic flow (i.e. where X{pe,k) > 0). In contrast to granular flow (cf. Fig. density 
waves develop at medium densities pe and small absolute wave numbers |fc| 7^ 0. At 
moderate densities, the absolute wave number related to the largest growth rate increases 
with density, corresponding to a decrease of the wave length of forming stop-and-go waves. 
As expected, traffic flow is only stable at low densities (free flow) and extreme densities 
(slow-moving traffic). The instabilities for large absolute wave numbers \k\ > 130/km 
would be connected with stop-and-go waves that move in forward direction with group 
velocity c{pe,k) > relative to Ve{pe) (b). However, these instabilities are physically 
irrelevant, since they lie in front of the lines ±27rpe which characterize the maximum wave 
numbers that can be propagated by vehicles with an average distance of 1/pe (= minimum 
wave length). The relevant instabilities at small absolute wave numbers \k\ < 2iTpc are 
connected with stop-and-go waves that move in backward direction with group velocity 
|c(pe5 k)\ > relative to K(pe)- In contrast to Fig. ^, the propagation speed has the right 
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